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Abstract 

The Jupiter Icy Moons Orbiter (JIMO) Space System is part of the NASA's Prometheus Program. As 
part of the JIMO engineering team at NASA Glenn Research Center, the structural design of the JIMO 
Heat Rejection Subsystem (HRS) is evaluated. An initial goal of this study was to perform sensitivity 
analyses to determine the relative importance of the input variables on the structural responses of the 
radiator panel. The desire was to let the sensitivity analysis information identify the important parameters. 
The probabilistic analysis methods illustrated here support this objective. 

The probabilistic structural performance evaluation of a HRS radiator sandwich panel was performed. 
The radiator panel structural performance was assessed in the presence of uncertainties in the loading, 
fabrication process variables, and material properties. The stress and displacement contours of the 
deterministic structural analysis at mean probability was performed and results presented. It is followed 
by a probabilistic evaluation to determine the effect of the primitive variables on the radiator panel 
structural performance. Based on uncertainties in material properties, structural geometry and loading, the 
results of the displacement and stress analysis are used as an input file for the probabilistic analysis of the 
panel. 

The sensitivity of the structural responses, such as maximum displacement and maximum tensile and 
compressive stresses of the facesheet in x and y directions and maximum VonMises stresses of the tube, 
to the loading and design variables is determined under the boundary condition where all edges of the 
radiator panel are pinned. Based on this study, design critical material and geometric parameters of the 
considered sandwich panel are identified. 


Introduction 

The structural design of the spacecraft has three phases: 

(1) Conceptual Design: 

(a) The feasibility studies are conducted. Design details related to geometry and material selection 
are identified. 

(b) Cost and risk is estimated 

(c) Requirements are derived 

(d) Candidate types of structures, materials are identified. 

(2) Preliminary design 

(a) The design parameters of the preferred candidates are recognized. 

(b) Material types and forms are selected. 

(c) Manufacturing plan is started. 

(d) Testing plans are started to be developed. 
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(3) Detail Design 

(a) The detailed final design properties are defined and final drawings for manufacturing are 
released. 

(b) Manufacturing and testing plans are developed. 

The intention of our study was to support the feasibility assessment of the JIMO Space System design. 

The HRS accepts the waste heat from the Power Conversion Subsystem (PCS) and transports the heat 
to radiator panels where it is rejected to space. Details about PCS for JIMO are given in references 1 
and 2. Therefore, the structural design of the HRS radiator panel is derived by the thermal design 
requirements to reject the required heat load produced in the spacecraft and maintain acceptable operating 
temperatures. Heat sources on the spacecraft come from environmental inputs and internal heat 
generation. The produced amount of heat is rejected to the space through HRS radiator panels. The 
radiator structural requirements should not exceed the available spacecraft heat rejection surface area. The 
design parameters of the radiator panel such as required diameter of the tubing depend upon the heat 
rejection capability of the radiator panel and the configuration of the tubing on the panel. Therefore 
design parameters are subject to change until the design uncertainties have been addressed. That is why it 
is important to asses what design parameters have significant affect on the structural behavior of the 
radiator panel under the design or manufacturing uncertainties in order to achieve more reliable and 
robust design. 

Most of the time analysts and designers are confronted with numerous uncertainties and design 
variability. To address these variables in traditional (deterministic) design, we apply safety factors. A 
reliability based design approach using probabilistic methods can explicitly model these uncertainties as 
random variables. When there is an uncertainty in the design or material parameters, using statistical 
techniques instead of comparing single numerical or test results, will help us to compare sets of data that 
reflect the variability of the structural responses as the result of uncertainties. 

For complicated and expensive aerospace structures, the ability to perform reliability and risk 
assessment without using extensive hardware testing is critical to design and certification. The use of 
computational simulation is relied upon increasingly more as performance requirements for engineered 
structures increase and as a means of reducing testing. Since structural performance is directly affected by 
uncertainties associated with models or in physical parameters and loadings, the development and 
application of probabilistic analysis methods suitable for use with complex numerical models is needed. 
The traditional established method of accounting for uncertainties is based on a factor-of-safety approach. 
However, this approach generally leads to over-conservative and hence non-optimized. 

The use of probabilistic methods will guide us to answer the following questions: 

(1) How these design uncertainties or manufacturing uncertainties affect structural performance of 
the HRS radiator panel such as maximum displacement, maximum tensile and compressive stresses? 

(2) Which input parameters will have significant affect on the structural responses? 

(3) Which input parameters should be addressed to achieve a reliable design and improve the 
quality? 


Geometry 

The outer dimensions of the radiator panel considered here are based on the preliminary radiator 
design feasibility studies. The radiator panel dimensions are shown in figure 1, and are 0.7087 in. 

(1.8 cm) panel thickness, 8 in. (20.32 cm) of panel width (direction normal to tube direction) and 6 in. 
(15.24 cm) of panel depth. The tube and saddle are located at the center of the panel. The tube rests in the 
saddle and the saddle is placed between two facesheets. The thickness of the facesheets is 0.025 in. 

(0.127 cm). The material used for facesheets is K1 100 Carbon-Carbon composite. The tube material is 
titanium, the outer diameter of the tube is 0.37 in. (0.9652 cm) and the thickness of the tube is 0.004 in. 
(0.01016 cm). 
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Material Properties 

The material and mass properties of the HRS radiator sandwich panel for GRC design are given in 
Table 1. Facesheets of the sandwich panel are made out of K1 100 Carbon-Carbon (C-C) composites and 
titanium is used for the heat pipe. The saddle and filler is constructed from Poco Foam and Reticulated 
Vitreous Carbon (RVC), respectively. 


Finite Element Model 

The finite element model of the HRS radiator sandwich panel was generated using MSC/Patran. The 
finite element model of the 3D sandwich panel as shown in figure 2, has 4,080 shell elements and 5,840 
hexagonal solid elements. The total number of nodes and elements are 18786 and 19920, respectively. 

The total weight of the 3D full sandwich panel section is 0.263 lb (0.12 kg). The sandwich panel was 
subjected to a 20 g loading and all sides of the panel were placed in a pinned boundary condition. 

Finite Element Analysis (FEA) was conducted using ANSYS version 8.1(ANSYS, Inc., Canonsburg, 
PA). The displacement contours of the sandwich panel, stress contours of the C-C facesheets in x and y 
directions and VonMises stress contours of the tube at mean probability are shown in figures 3 to 6. 

The Probabilistic Model 

This study shows how FEA and probabilistic design can be used to simulate the effects of design or 
manufacturing uncertainties on the structural responses such as maximum displacement, maximum tensile 
and compressive stresses. The random variables describe modeling and manufacturing uncertainty. The 
uncertainties may come from conceptual modeling because of lack of data on the physical process 
involved or lack of system knowledge or from material properties such as elastic properties, density, local 
imperfections, etc. A Monte Carlo method is used for probabilistic simulation and Latin Hypercube 
method used for sampling (refs. 3 and 4). 

The Monte Carlo Simulation method lets you simulate how virtual components behave the way they 
are built. Each one of the simulation loops represents one manufactured component that is subjected to a 
particular set of design or load conditions. The Latin Hypercube Sampling (LHS) technique has a sample 
“memory,” meaning it avoids repeating samples that have been evaluated before (it avoids clustering 
samples). It also forces the tails of a distribution to participate in the sampling process. For the Latin 
Hypercube Sampling technique the range of all random input variables is divided into n intervals with 
equal probability, where n is the number of sampling points (ref. 4). For each random variable each 
interval is “hit” only once with a sampling point. The process of generating sampling points with Latin 
Hypercube restricts the sampling points within the respective interval. 

For this particular application, 100 finite element simulations using randomly generated input values 
provided sufficient data to accurately estimate the cumulative distribution function (CDF). The 
probabilistic analysis method repeatedly evaluates the deterministic model to generate samples of the 
model response, from which the response statistics and probability are approximated. The probabilistic 
analysis process in ANSYS version 8.1 consists of the general steps that are shown in figure 7 (ref. 3). 

Table 1 contains a description of the 7 random variables. All variables were allowed to vary using 
normal distribution. The scatter range of the uncertainties considered for the probabilistic evaluation is 
defined as percentage or coefficient of variation from the mean value. The standard deviations for all of 
the parameters are established based on engineering judgment. The type of distribution considered for all 
parameters is Normal distribution. The mean value and standard deviation of all random input variables 
are shown in Table 2. 
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The Probabilistic Analysis Results 

The effect of the uncertainties in the design variables on the maximum displacement, maximum 
longitudinal and transverse tensile stresses of the C-C composite facesheets and on the minimum 
longitudinal and transverse compressive stresses of the C-C composite facesheets and on the maximum 
VonMises equivalent stress of the heat pipe is evaluated. The Cumulative Distribution Function (CDF) 
and sensitivity factors are determined. 

The plots of the CDFs for the maximum displacement, maximum longitudinal tensile and 
compressive stress, maximum transverse tensile and compressive stress and maximum VonMises stress 
are shown in figures 8 to 13. The upper and lower lines shown in CDF plots display the confidence 
bounds. The 95 percent confidence level is used for all responses. The confidence bounds quantify the 
accuracy of the probability level. 

The CDF for the maximum displacement is shown in figure 8. The mean probabilistic maximum 
displacement is 0.048 in. The 0.01 percent probabilistic displacement is 0.075 in. and 99.99 percent 
probabilistic displacement is 0.025 in. The CDF for the longitudinal tensile and compressive stresses in 
the C-C composite facesheet are shown in figures 9 and 10. The mean maximum longitudinal tensile 
stress in the C-C composite facesheet is about 14677 psi with the 0.01 percent probabilistic stress of 6000 
psi and 99.99 percent probabilistic stress of 27900 psi. The mean probabilistic maximum longitudinal 
compressive stress of the C-C composite facesheet is about -14677 psi with the 0.01 percent probabilistic 
stress of -27900 psi and 99.99 percent probabilistic stress of -6000 psi. The scatter in the maximum 
longitudinal tensile or compressive stress is about 21900 psi. 

The CDF for the transverse tensile and compressive stresses in the C-C composite facesheet are 
shown in figures 1 1 and 12. The mean probabilistic maximum transverse tensile stress in the C-C 
composite facesheet is about 18775 psi with the 0.01 percent probabilistic stress of 12500 psi and 
99.99 percent probabilistic stress of 26300 psi. The mean probabilistic maximum transverse compressive 
stress of the C-C composite facesheet is about -18775 psi with the 0.01 percent probabilistic stress of 
-26300 psi and 99.99 percent probabilistic stress of-12500 psi. The scatter in the maximum transverse 
tensile or compressive stress is about 13800 psi. 

The CDF for the maximum VonMises stress in the heat pipe shown in figure 13. The mean 
probabilistic maximum VonMises stress is about 149770 psi with the 0.01 percent probabilistic stress of 
74000 psi and 99.99 percent probabilistic stress of 224000 psi. The scatter in the maximum VonMises 
stress is about 21900 psi. 

The probabilistic sensitivities of the input parameters on the maximum displacement of the sandwich 
panel is shown in figure 14. The most important input parameters that have significant effect on the 
maximum displacement, because of the uncertainties, are applied g load, facesheet thickness, facesheet 
density and facesheet Poisson’s ratio, in importance order. The maximum longitudinal tensile and 
compressive stress of the C-C composite facesheet is most sensitive to uncertainties in the facesheet 
thickness, g load, longitudinal modulus of elasticity of the C-C composite facesheet, facesheet density and 
in-plane Poisson’s ratio of the C-C composite facesheet, respectively (figs. 15 and 16). The uncertainties 
in g load and facesheet density have significant affect on the transverse maximum tensile and 
compressive stresses of the C-C composite facesheet (figs. 17 and 18). The maximum VonMises stress of 
the heat pipe is most sensitive to uncertainties in the g load, facesheet thickness, facesheet and density 
(fig. 19). 

Positive sensitivity values indicate a positive change in the mean value will result in an increase in the 
computed probability. Negative sensitivities indicate that a positive increase in the mean value will result 
in a decrease in the computed probability. 
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Conclusions 


For the JIMO HRS radiator panel section, the probabilistic evaluation of the structural responses to 
the input parameters is presented. The Cumulative Distribution Function (CDF) and response sensitivity 
factors are calculated by the computational simulation in ANSYS version 8.1. The effect of the 
uncertainties in the input parameters on the maximum displacement, maximum longitudinal and 
transverse stresses of the facesheet and maximum VonMises stress of the heat pipe are evaluated. The 
following conclusions are drawn: 

1. Controlling or modifying the scatter in the applied g load, facesheet thickness, facesheet density 
and in-plane Poisson’s ratio of the C-C composite facesheet is critical to control the maximum deflection 
of the sandwich panel. 

2. In general, results indicated that uncertainties of the facesheet thickness, g load, longitudinal 
modulus of elasticity of the C-C composite facesheet, facesheet density and in-plane Poisson’s ratio of the 
C-C composite facesheet have a significant effect on the maximum longitudinal tensile and compressive 
stresses of the facesheet. 

3. The uncertainties in g load and facesheet density dominate and have significant effect on the 
maximum transverse tensile and compressive stresses of the facesheet. 

4. The maximum VonMises stresses of the Heat Pipe is sensitive to the uncertainties in the g load, 
facesheet thickness and facesheet density. 

5. The computational capability makes probabilistic analysis acceptable for practical applications. 

6. Probabilistic sensitivities can be used to provide guidance for necessary design changes. 

7. Probabilistic methods can be used to quantify the reliability of the HRS radiator panel to achieve a 
more robust design and improve the quality. 
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TABLE L— MATERIAL PROPERTIES OF THE GRC RADIATOR SANDWICH PANEL 


Properties 

Carbon/carbon 

(orthotropic) 

Titanium 

(isotropic) 

POCO foam 
(isotropic) 

RVC foam 
(isotropic) 

Eli, psi 

6.3091xl0 7 

1.52xl0 7 

5.8 xio 4 

6xl0 3 

E 22 , psi 

1.6824xl0 7 

n/a 

n/a 

n/a 

E 33 , psi 

1.3050xl0 7 

n/a 

n/a 

n/a 

V12 

0.21 

0.31 

0.4082 

n/a 

V23 

0.04 

n/a 

n/a 

n/a 

V31 

0.04 

n/a 

n/a 

n/a 

G 12 , psi 

2.84x1 0 6 

n/a 

n/a 

4.4x1 0 3 

G 23 , psi 

9xl0 5 

n/a 

n/a 

n/a 

G 31 , psi 

9xl0 5 

n/a 

n/a 

n/a 

p, lb/in. 3 

0.0692 

0.163 

0.018449 

0.001794 


TABLE 2.— INPUT RANDOM VARIABLES FOR RADIATOR SANDWICH PANEL 


Parameters 

Distribution 

Mean 

Standard deviation, 
percent 

Units 

El 

Normal 

6.3091x107 

12 

psi 

E2 

Normal 

1.6824x107 

12 

psi 

Nul2 

Normal 

0.21 

12 


FSTHICK 

Normal 

0.025 

12 

in. 

TUBETHICK 

Normal 

0.0039 

12 

in. 

DENSITY 

Normal 

0.0692 

12 

lb/in. 3 

GLOAD 

Normal 

20 

12 

g 



3.671 in. 

0.658 in. 

3.671 in. 

^ w 




l 


: 0.025 in. 


T 


0.6587 in. 


10.025 in. 


Figure 1. — Cross section of the 8 by 6 in. sandwich panel. 
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Figure 2. — The top and side views of the finite element model of the HRS radiator sandwich panel. 


NODAL SOLUTION 


STEP=1 
SUE =1 
TIME=1 

UZ (AVG) 

RSYS=0 

DMX = . 048472 
SHN =-.043472 


AN 



Figure 3. — Displacement contours of the HRS radiator sandwich panel at mean probability. 
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NODAL SOLUTION 


AN fS 



SUE =1 
TIHE=1 

SX (AVG) 

RSYS=0 

DHX = . 04847Z 
SUN =-15879 
SHX =15879 


Figure 4. — Longitudinal (in x direction) stress contours of the C-C composite 
facesheets at mean probability. 



Figure 5. — Transverse (in y direction) stress contours of the C-C composite 
facesheets at mean probability. 


NASA/TM— 2006-2 14 1 1 6 



Figure 6. — VonMises stress contours of the titanium heat pipe at mean probability. 
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Choose the probabilistic design tool or method 


Execute the loops required for the probabilistic design analysis 


Review and plot the results of the probabilistic analysis 

Figure 7. — The probabilistic analysis process in ANSYS Version 8.1. 
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9 9 . 99 



AN 
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SKEW -0.42457E+00 
KURT 0 . 55103E+00 
MIN -0.73875E-01 
MAX -0.27262E-01 


Conf idsric* Li in it 
95 . 00 % 


Figure 8. — CDF of the maximum displacement (in inches). 
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Result Set RAD IATOR- P AWE L 


99 . 99 



ANSYS 

HE AN 0.14577E+05 
STDEU 1 0 . 36873E+04 
SKEW 0.77845E+00 
KURT 0.53499E+00 
MIN 0.50378E+04 
MAX 0 . 27590E+05 

CorL£id*rict Limit 
95.00% 


Figure 9. — CDF of the maximum longitudinal tensile stress (in psi) acting on the C-C composite facesheet. 
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ANSYS 


Cumulative Distribution Function 
Result Set RAD IATOR- P AWE L 


9 9 . 99 



HEAN -0.14577E+O5 
3TDEU 0.36874E+04 
3KETO -0.77845E+O0 
KURT 0.53499E+00 
HIN -0.27590E +05 
MAX -0.5037SE+04 

Conf idsric* Li in it 
95 . 00 % 


Figure 10. — CDF of the maximum longitudinal compressive stress (in psi) 
acting on the C-C composite facesheet. 


Cumulative Distribution Function 
Result Set RAD IATOR- P AWE L 


99 . 99 



AN 


HEAN O.18775E+05 
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CorL£id*rict Limit 
95.00% 


Figure 11. — CDF of the maximum transverse tensile stress (in psi) 
acting on the C-C composite facesheet. 
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Cumulative Distribution Function 
Result Set RAD IATOR- P AWE L 


9 9 . 99 



ANSYS 


KEAN -0.13775E+05 
3TDEU 0.26193E+04 
3KETO -0 . 20159E+00 
KURT 0.44433E+00 
KIN -0.26515E +05 
HAX -0.12509E+05 


Conf idsric* Li in it 
95 . 00 % 


Figure 12. — CDF of the maximum transverse compressive stress (in psi) 
acting on the C-C composite facesheet. 


Cumulative Distribution Function 
Result Set RAD IATOR- P AWE L 


99 . 99 



AN 

HE AN 0.14977E+06 
3TDEU 0.25729E+05 
3KETO 0.50855E+00 
KURT 0.71342E+00 
WIN 0 . 79950E+05 
MAX 0 . 22331E+06 

CorL£id*rict Limit 
95.00% 


Figure 13. — CDF of the VonMises equivalent stress (in psi) acting on the heat pipe. 
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Figure 14. — Sensitivities of the maximum displacement. 
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Figure 15. — Sensitivities of the maximum longitudinal tensile stress of the C-C composite facesheet. 
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Figure 16. — Sensitivities of the maximum longitudinal compressive stress of the C-C composite facesheet. 
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Figure 17. — Sensitivities of the maximum transverse tensile stress of the C-C composite facesheet. 
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Figure 18. — Sensitivities of the maximum transverse compressive stress of the C-C composite facesheet. 
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Figure 19. — Sensitivities of the maximum VonMises stress of the heat pipe. 
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